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BiFeC>3 thin films epitaxially grown on SrRuOg-buffered (OOl)-oriented SrTiC>3 substrates show 
orthogonal bundles of twin domains, each of which contains parallel and periodic 71° domain walls. 
A smaller amount of 109° domain walls are also present at the boundaries between two adjacent 
bundles. All as-grown twin walls display enhanced conductivity with respect to the domains during 
local probe measurements, due to the selective lowering of the Schottky barrier between the film 
and the AFM tip (see S. Farokhipoor and B. Noheda, Phys. Rev. Lett. 107, 127601 (2011)). In this 
paper we further discuss these results and show why other conduction mechanisms are discarded. In 
addition we show the crucial role that oxygen vacancies play in determining the amount of conduction 
at the walls. This prompts us to propose that the oxygen vacancies migrating to the walls locally 
lower the Schottky barrier. This mechanism would then be less efficient in non-ferroelastic domain 
walls where one expects no strain gradients around the walls and thus (assuming that walls are not 
charged) no driving force for accumulation of defects. 

PACS numbers: 77.80.Dj, 68.37.Ps, 73.30.+y, 73.61.Lc 



I. INTRODUCTION 

The importance of domain walls and their influence 
in the response of materials have been recognized since 
decades p} but it has been mainly magnetic domains 
walls, which typically extend to decades of unit cells, the 
ones most thoroughly investigated. Less widely known 
was, perhaps, that ferroelastic twin walls, spanning a few 
unit cells in width, are exceptional physical nano-objects 
with remarkable properties [2, 3|, being a notable example 
the observation of superconductivity at the ferroelastic 
domain walls of doped WOajij]. Indeed, ferroelastic, or 
twin, walls are also inherent symmetry-breaking entities 
within a given material, with the potential to display the 
distinct additional physical responses associated to that 
fact. Moreover, they also represent areas of intense strain 
gradients that can attract large concentrations of dopants 
or defects, which will also modify greatly and very lo- 
cally the materials' properties [f|. In addition, the strain 
gradients associated to those walls can induce electrical 
dipoles, due to the so-called flexoelectric effect |6r9]- 

Even if tantalizing results on twin walls were known 
and twin walls had been accessed and imaged [13, [HI > 
it is only recently that it is made clear that, using thin 
film deposition techniques, reproducible control and tun- 
ability of the type, orientation and periodicity of do- 
main walls, separated by distances of only tens to hun- 
dreds of nanometers, is possible. Thus, also due to the 
widespread use of local probes, it is now realistic to think 
of usin g th e domain walls and their unique properties as 

devices|H,[ll. 

But what has boosted the general interest on do- 
main walls is the recent discovery that domain walls of 
multiferroic (ferroelectric and antiferromagnet) BiFeC>3 
thin films are considerably more conducting than the 
domains and, thus, provide well-defined, local paths of 
conduction through the thin films[14]. BiFeC>3 thin 



films are monoclinic[15[ but, in fact, they are quasi- 
rhombohedrally distorted and, thus, the ferroelastic do- 
main walls that they present are very similar to those ex- 
pected in a rhombohedral perovskite, that is 71° and 109° 
domain walls [16]. Initially, conducting behavior was ob- 
served only in artificially- written 109° domain walls 14]. 
However, it was later reported that also as-grown 71° 
walls displayed enhanced conductivity with respect to 
the domains 17J , although that was one order of mag- 
nitude smaller than the conductivity found in as-grown 
109° domain walls. This seems to be in agreement with 
the theoretical calculations that assign the origin of con- 
duction to a decrease in the band gap of the material 
at the domain walls [3, [l7j]- However, in samples grown 
under different conditions, we have observed that con- 
ductivity at 71° domain walls can be as large as that in 
109° domain walls and that its main origin is the lowering 
of the Schottky barrier between the BiFeC>3 n-t ype semi- 
conducting film and the metallic top electrode |18j. The 
predicted reduction of the band gap at the walls would 
then appear as a secondary effect, responsible for subtler 
changes in conduction. 

Therefore, if we want to achieve control of walls ei- 
ther as devices or as parts of a device, it is important 
to investigate the mechanism or mechanisms for selec- 
tive conduction through domain walls and, in case that 
there is more that one possibility, to find out which one 
is most robust and convenient. This becomes now essen- 
tial because, according to recent reports, conduction at 
domain walls is a rather general phenomenon[19], only 
the most relevant mechanisms seem to vary. As one 
can expect, different conduction mechanisms are likely 
to have different relevance in different types of walls. 
Fowler-Nordheim (FN) tunneling has been put forward 
as the mechanism for conduction in artificially-written 
109° domain walls [20l [2l| . Recent work has also focused 
on walls between ferroelectric 180° domains, which do 
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not involve strain gradients [lj|[22[ and indeed a different 
conduction mechanism (bulk-limited Poole-Frenkel emis- 
sion) has been observed in this case. Similarly, charged 
walls show again different conduction characteristics |23| . 
Later work has gone further in exploring the conduction 
involved in truly one-dimensional conduction paths de- 
fined by vortex states in ferroelectrics, around which both 
charged walls and large strain gradients coexist [24] 

In this paper we focus on ferroelastic domain walls of 
BiFeC>3 thin films to complement the results reported in 
ref.[l8[. We include here a more detailed description of 
our samples and a more detailed analysis of the differ- 
ent conduction mechanisms that have been considered 
and the reasons why they have been discarded, as well 
as how the Schottky barriers have been calculated. We 
have also pointed out[IH that the lowering of the bar- 
rier at the walls is likely to be induced by migration of 
charged defects to the walls (as it is common in other 
perovskites[5|), since that would create a potential step 
at the interface 25] . We postulated that the charged de- 
fects are oxygen vacancies, based on the observed en- 
hancement of current in samples with increased oxygen 
vacancy concentration. In this paper we show evidence 
of the crucial role played by the oxygen content in deter- 
mining the domain wall conductivity. 



II. EXPERIMENTAL 

BiFe03 thin films were grown by pulsed laser deposi- 
tion (PLD) assisted by RHEED using a system designed 
by Twente Solid State Technology (TSST)|6| provided 
with an excimer laser (Lambda Physik COMPex Pro 
205) filled with KrF to produce a wavelength of 248 nm. 
The operating laser frequency and fluence were 0.5Hz 
and 2 J/cm 2 , respectively. The chamber was evacu- 
ated up to a background pressure of 10 -8 mbar. The 
films were grown at a temperature of 670°C in an oxy- 
gen pressure of 10 _1 mbar. After growth, the oxygen 
pressure was increased to lOOmbar and the films were 
cooled down at a rate of 3° C/min, except in particu- 
lar cases in which we intended to investigate the effect 
of different oxygen vacancy content, for which the cool- 
ing rate was decreased(increased) or the oxygen pressure 
was increased(decreased). The films are grown on single- 
terminated (TiO2)(001)-oriented SrTiC>3 substrates [27j 
covered by a buffer electrode layer of SrRu03. The 5nm 
thick SrRuC>3 layer was deposited also by PLD imme- 
diately before the BiFe03 deposition, with a substrate 
temperature of 600°C, an oxygen pressure of 0.13 mbar 
and laser frequency and fluence of 1 Hz and 2 J/cm 2 , re- 
spectively. 

(OOl)-oriented BiFeC>3 films with thickness ranging be- 
tween 40 and 68 nm are used in this study. These films 
show atomically fiat terraces (see Figure la-b, flat inter- 
faces and good crystal quality (Figure lc). They form 
monoclinic, pseudo-rhombohedral, domains with the po- 
larization very close to the pseudo-cubic [lll] p direction. 
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FIG. 1: a) 2/j.m x 2^tm atomic force microscope (AFM) image 
of the topology of a 50nm thick BiFeC>3 film grown by PLD 
on a 5nm thick SrRuOa electrode on a low miscut, single- 
terminated, (OOl)-oriented, SrTiOs substrate, b) A linear 
scan of such image shows atomically flat terraces following 
the substrate steps, c) XRD pattern of the same film around 
the (001) reflection of the BiFeOa film (low-angle peak) and 
SrTiC>3 substrate (high-angle peak) The thickness fringes cor- 
responding to both the film and the electrode are clearly visi- 
ble and indicate the excellent quality of the interfaces. In the 
inset, a reciprocal space map around the pseudo-cubic (113) 
reflection of the film shows that it is monoclinic with the [001] 
direction parallel to that of the substrate [31]]. 

From the eight possible polarization orientations, only 
the four down oriented domains are present in these films, 
in agreement with other reports[28, 29], and thus no 180° 
domain walls are present. The observed domains are or- 
ganized in bundles [30]] of parallel 71° domain walls [3lj]. 
Adjacent bundles are oriented orthogonal to each other 
and at the boundaries in between bundles some 109° do- 
main walls are also found. 

From the reciprocal space maps around non-specular 
reflections (inset of Figure lc and ref.[3l|), it can be 
inferred that the (001) atomic planes of the film are 
parallel to those of the substrate and buffer layer and, 
thus, that there is no out-of-plane tilt or buckled inter- 
face, as reported for some (OOl)-oriented rhombohedral 
perovskites [28i [32l .331] . This is important because the na- 
ture of the interfaces will largely influence the conduction 
properties of the films [341] 

We have performed DC transport measurements on a 
VEECO (now Bruker) Dimension V Conductive Atomic 
Force Microscope (CAFM) provided with a (Co-Cr 
coated Si) metallic tip. Electrically, the tip is connected 
to the ground of the microscope and a bias voltage is ap- 
plied to the metallic sample holder, which is connected 
to the SrRuC>3 bottom electrode using silver paste on the 
side of the sample. A so-called TUNA™ amplifier al- 
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lows to measure currents through the film in the range of 
100 fA-10/xA, using four different gain factors (and sensi- 
tivities) . Measurements are performed both as a function 
of temperature and applied voltage in the range of -5 V 
to +5V, well below the voltages at which currents asso- 
ciated with domain switching are observed (at about 7 
V). The temperature was increased in steps of 10 K in 
ambient conditions from room temperature up to 155 °C. 
In order to prevent spurious transient currents that are 
related to ionic conductivity, we included a delay time 
of several minutes from the application of the voltage to 
the starting of the measurements. 



III. RESULTS AND DISCUSSION 

Figure 2c shows localized areas of enhanced conduction 
around both 71° and 109° domain walls in the as-grown 
(001)-BiFeO3 films described above, as we have previ- 
ously reporte dll8ll . Similar results have been observed 
by Chiu et al.|17j using scanning tunneling microscopy. 
However, in our films, unlike in those of ref.|17j|. the mag- 
nitude of the current in the 71° walls is as large as that of 
the 109° walls. This difference between the two studies 
could be due to the nature of the interfaces (out-of-plane 
twinning versus in-plane twinning) [32l l34T |. which may 
invoke different conduction mechanisms, as mentioned in 
the previous section. Moreover, the 71° domain walls in 
our investigated films showed a somewhat more (statis- 
tically) consistent behavior than that of the 109° walls, 
probably due to the fact that the later are located at the 
edges of the orthogonal bundled regions, where the strain 
is enhanced, possibly inducing dislocations or other de- 
fects. Therefore, we have focussed most of our studies on 
the 71° domain walls. 

The magnitude of the measured current at the walls 
is of about 10 pA at 4 V at room temperature and can 
increase up to 3nA for temperatures around 150°C. By 
measuring I-V curves as a function of temperature (above 
room temperature), we have evaluated the three most 
likely mechanisms leading to conduction through a metal- 
semiconductor-metal system, as well as the possibility of 
space charge limiting the current. We have considered a 
linear dependence between the applied voltage and the 
electric field, which is not obvious giving that one of the 
electrodes is a sharp tip. However, recently, Guyonnet et 
al.jlS] have shown that this is, indeed, a good approxi- 
mation. 

Figure 3a) represents the data in coordinates related 
to Fowler-Nordheim tunneling, that is ln(I/V 2 ) is plotted 
against 1/V for different temperatures. Fowler-Nordheim 
tunneling should be temperature-independent and, in 
these coordinates, should show a linear dependence. Fig- 
ure 3a shows quite linear behavior mainly in the 3-5V 
region but there is a clear temperature dependence that 
indicates that this is not the relevant mechanism for con- 
duction. 

Another possible mechanism is Poole-Frenkel (PF) 
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FIG. 2: Piezo-force microscopy (PFM) amplitude (a) and 
phase (b) images of a BiFe03 film as those described in the 
text. Both images are needed to determine the type of domain 
walls, c) Conductive-AFM (C-AFM) image of the same area 
of the film showing enhanced currents at the domain walls. 
Both 71° and 109° domain walls show similar magnitudes of 
the current. 
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FIG. 3: The measured I-V curves for different temperatures at 71° domain walls are plotted with different axes in order to 
make the curves linear and temperature-independent for a) Fowler-Nordheim (FN) tunneling; b) Poole-Frenkel (PF) emission; 
c) Simmons-Richardson-Schottky (SRS) emission and d) space charge limited (SCL) current, a) shows temperature dependent 
slopes that are inconsistent with the FN mechanism; the slope in b) leads to a too small optical dielectric constant of 0.4±0.1; 
while the slope of c) gives the right value of 6.5±1.5. d) shows the data for voltages below 3.7 V, which deviate from linearity 
in figure 3c). Power laws with exponents close to 2 seem to fit the data reasonable well 18] in agreement with a space charge 
limited regime. 



emission, by which current is supported by thermal emis- 
sion from trapped carriers in the bulk of the semicon- 
ductor, where the trap barrier is lowered by the electric 
field. This has been reported to be the main conduc- 
tion mechanism in some ferroelectric oxides, including 
BiFe0 3 [H and PZtJT^. The field and temperature de- 
pendence of the PF currents are similar to those of the 
Richardson-Schottky (RS) emission, a third mechanism 
to be considered. The best way to distinguish between 
the two of them is by extracting the high-frequency di- 
electric constant from the slopes of the linearized plots 
and see in which case it agrees with that known for 
the material. [Hj] . In addition, the Richardson-Schottky 
equation is only strictly correct in the case of semicon- 
ductors with relatively small electronic mean free paths; 
while for insulators the modification proposed by Sim- 
mons (which we will call Simmons-Richardson-Schottky 



(SRS) emission) [361] should apply. 

Thus, in Fig 3 b)-c), we plot the data in the coor- 
dinates that make the current data due to PF emis- 
sion or SRS emission linear and temperature- independent 
(Tln(I/V) versus V 1 / 2 and Tln[I/(VT 3 / 2 )] versus V 1 / 2 , 
respectively). In this way, in both plots, all the curves 
(independent of temperature) should show the same in- 
tercept and the same slop e, which is inversely-related to 
the dielectric constant j35|. The linearity of the fits is con- 
siderably improved compared to the FN tunneling fit, in 
particular at higher temperatures and higher voltages. 
One can see that the slopes are in both cases (PF as 
well as SRS) temperature-independent in a very good 
approximation. The collapse of the intercept seems bet- 
ter in the SRS case (leading to a well-defined Schottky 
barrier). However, it is the value of the optical dielectric 
susceptibility extracted from the slopes what makes us 
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FIG. 4: Ln(I/(VT 3/2 ) versus T _1 for two particular voltages, 
V= 3.8 V (a) and V= 4.8V (b), limits within which the SRS 
mechanism fulfills. The slopes of those curves give rise to what 
has been called (j> avV arent [3J]. c) According to the SRS mech- 
anism, 4>apparent vs. V 1//2 should show linear behavior with 
the intercept at the origin being the height of the Schottky 
barrier (<f>o)- The <f> values extracted in this way for domains 
and domain walls are indicated in the plot. 



conclude that Schottky emission (SRS) is the pertinent 
conduction mechanism in this case: In the case of PF 
emission a value of 0.4±0.1 for the optical dielectric con- 
stant is found from the slopes of Fig 3b) for 71° domain 
walls and a value of 0.16±0.06 is found for the domains 
(figure not shown). For SRS emission, a dielectric permit- 
tivity of 6.5±1.5 is found at the 71° domain walls (from 
figure 3c), in perfect agreement with the expected value 
of 6.25 (see e.g. ^3). A value of 3.6±0.1 is found at 109° 
domain walls; while a worse agreement of 1.8±0.3 is ob- 
tained for the domains. It is worth to mention that, due 
to the low values of the current inside the domains, these 
data are very noisy and only a few higher temperature 
curves could be analyzed. 

To calculate the Schottky barrier height we follow 
ref. (H: we first plot Ln(I/VT 3 / 2 ) versus 1/T. A 
linear fit should be obtained for the proper voltage 
regime (voltage large enough to support the Schot- 
tky emission mechanism, in our case above 3.8V). Fig- 
ures 4 a)-b) show the fits for two specific (the largest 
and the lowest) voltages, ^apparent is defined as the 
slope of these linear fits. This is the exponential term 
in the Simmons- Richardson-Schottky emission (36|: (- 
$ /ki3) + (c/k i 3)(eV/d7re K) 1 / 2 , where $ Q is the Schot- 
tky barrier, kg is the Boltzmann constant, e D is the per- 
mittivity of free space, K is the dielectric constant, V is 
the applied DC bias and d is the thickness of the film. 
Next, we plot the slopes versus V 1 / 2 (Figure 4c). The 
intercept of the obtained graph (at V=0) is the Schottky 
barrier height. The barrier height for 71° domain walls is 




FIG. 5: C-AFM images on an 2/imx 2/im area for three sam- 
ples grown under identical conditions but cooled down after 
growth a) under an oxygen pressure of P02= 300 mbar and a 
cooling rate of 2K/min; b) under P02= lOOmbar and cooling 
rate of 3 K/min and c)under P02= lOOmbar and cooling rate 
of 30-4 K/min. The current scale is the same for the three 
images d) Linear scans across 1 /im of the previous images 
are shown in red (for a)), blue (for b) and black (for c). 



determined to be (0.8 ± 0.1) eV, however, inside domains 
this value changes to (2.6±0.7) eV. 

One can see in Fig. 3b-c) that at the lower voltages 
the data deviates from the linear trend indicating a differ- 
ent behavior. The lower the temperature, the higher the 
voltage at which the SRS regime crosses over to this low 
voltage regime. It is known that the general signature 
of space charge limited (SCL) conduction in an insula- 
tor of thickness d (with parallel-plate electrodes) is Ioc 
V n /d n+1 , where V is the applied voltage and n>2 (n= 2 
in the absence of traps) [39] . Typically, the SCL regime is 
found at the higher voltages (once there are enough car- 
riers in the insulator to build up space charge). However, 
charged defects or, in a ferroelectric material, the polar- 
ization surface charges themselves can behave similarly 
to the space-charge regions and could lead to a similar 
I-V dependence [401] . In Fig. 3d) we plot the low volt- 
age data that, according to Figure 3c), deviates from the 
linear Schottky emission behavior. As discussed before, 
this deviation is more frequently observed at the lower 
temperatures and it is in reasonable agreement with an 
n= 2 slope. In this low voltage regime (2.5-3.5V), Arrhe- 
nius behavior of the current has been observed, with a 
thermal activation energy of 0.73(9)eV[l8|]. This value is 
consistent with trapped electrons from clusters of oxygen 
vacancies in perovskites.[20l|. 

That oxygen vacancies play a very important role in 
controlling the conduction in these films has been clearly 
observed by making samples with different oxygen con- 
tents, cither by changing the cooling rate or the oxygen 
pressure after growth, but otherwise growing them un- 
der identical conditions. These results are summarized 
in Figure 5, where three C-AFM images are displayed for 
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samples with low, medium and high vacancy content. An 
increasing amount of current is obtained with decreasing 
oxygen content. A clear dependence of the current with 
growth oxygen pressure was also observed by Seidel et 
al.[H in 109° twin walls. 

All this brings forward the following scenario: the con- 
duction observed in the low voltage regime is coming from 
thermally excited electrons from defect states (possibly 
related to oxygen vacancies) located close to the bottom 
of the conduction band. With increasing voltages, oxy- 
gen vacancies move towards the surface and lower the 
Schottky barrier with the top electrode, which eventu- 
ally allows large conduction through the films. This ef- 
fect is, however, not homogeneous across the film: strain 
gradients associated to ferroelastic domain walls create 
an inhomogeneous distribution of oxygen vacancies and 
the selective reduction of the Schottky barrier around the 
walls. The mechanism for conduction is thus expected to 
be different in non-ferroelastic domain walls. 

To conclude, we have observed that in (001)- BiFcC>3 
thin films grown on SrRuOs-buffered (001)-SrTiO3 sub- 
strates, the currents above about 3.8 V (for films with 



thickness of about 50nm) are determined by the Schottky 
barrier between the BiFeOa film and the Co-Cr tip. This 
barrier is largely decreased in the domain walls with re- 
spect to the domain, causing the observed conduction en- 
hancement at the walls. Oxygen vacancies are observed 
to play a crucial role in this. At low voltages the current 
behaves like it is limited by space charge. The importance 
of the interplay between the surface polarization charges 
and the charged defects have been recently put forward 
in BiFe0 3 single crystals by[Il[ and in films [23ll24j. This 
interplay is, of course, very different at the domain and 
at the walls and it has been investigated in relation with 
domain/ domain wall conductivity in the present BiFeC>3 
thin films. The results will be discussed elsewhere [42} . 
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